Dynamics of periodic anticrossings: Decoherence, pointer states and hysteresis curves 
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We consider a strongly driven two-level (spin) system, with a periodic external field that induces 
a sequence of avoided level crossings. The spin system interacts with a bosonic reservoir which leads 
to decoherence. A Markovian dynamical equation is introduced without relying on the rotating 
wave approximation in the system-external field interaction. We show that the time evolution of 
the two-level system is directed towards an incoherent sum of periodic Floquet states regardless of 
the initial state and even the type of the coupling to the environment. Analyzing the time scale of 
approaching these time- dependent pointer states, information can be deduced concerning the nature 
and strength of the system-environment coupling. The inversion as a function of the external field 
is usually multi-valued, and the form of these hysteresis curves is qualitatively different for low 
and high temperatures. For moderate temperatures we found that the series of Landau-Zener- 
Stiickelberg-type transitions still can be used for state preparation, regardless of the decoherence 
rate. Possible applications include quantum information processing and molecular nanomagnets. 



I. INTRODUCTION 

Crossings and anticrossings of energy levels play im- 
portant role in various physical systems. The case when 
the level scheme is time dependent is of special inter- 
est, because the separation of the levels strongly influ- 
ences the dynamics, e.g., depending on the parameters, 
transitions may occur around an anticrossing. For a 
two-level system with a linearly time-dependent Hamilto- 
nian, the analytically solvable Landau-Zener-Stiickelberg 
(LZS) model [H, 0, S] reflects the most important proper- 
ties of the dynamics. Phenomena that can be described 
using this model include the dynamics of molecular vi- 
brations 0, multiphoton transitions and ionization 
of Rydberg atoms [1, 0, H, @, molecular nanomagnets 
[lOl [ill . [T3 ITst and also quantum information process- 
ing with superconducting qubits [H, [l^ . Multilevel and 
nonlinear generalizations of the LZS model have also 
been studied extensively (see e.g. [H, [13, [H, Ell ) . The 
LZS model itself can describe the adiabatic limit, when 
the system follows the instantaneous eigenstates of the 
Hamiltonian, as well as the case when a sudden transition 
takes place. Generally (between the two extreme situa- 
tions mentioned above) , an LZS transition splits the state 
of the system into two parts, which are almost orthogonal 
when the transition region is left. In this sense, consid- 
ering periodic driving that forces the system to return 
to the transition region after a half-cycle, fundamental 
which- way interference effects [13] can appear. This kind 
of interference is highly sensitive to dephasing and other 
decoherence mechanisms [2l|, which are the phenomena 
to be investigated in the present paper. 

We consider a two-level system driven by an external 
field leading to periodic anticrossings. Mathematically, 
this is similar to the LZS model with harmonic terms 
replacing the usual linear ones in the Hamiltonian. The 
spin 1/2 system is assumed to be embedded in a thermal 



reservoir, and using Floquet theory [24I we introduce the 
appropriate master equation in a systematic way. Inves- 
tigating the dynamics induced by the master equation 
we determine the direction and time scale of the deco- 
herence. We find that the time-dependent Floquet states 
play the role of pointer states [2^ . The temperature de- 
pendence of the process is also investigated with a special 
focus being on the population difference (that is, (cz)) as 
a function of the external field. In the context of molec- 
ular nanomagnets, similar magnetization curves have re- 
cently been investigated in pulsed fields both experimen- 
tally (see e.g. [Ulil,!!!]) and theoretically [2i,|27|. Our 
findings in the periodic case generalize the result that 
the multi-valuedness of the hysteresis curves in these sys- 
tems reflects the relaxation towards the quasistationary 
(equilibrium) solutions, which are related to single- valued 
magnetization curves. 



II. DYNAMICAL EQUATIONS 

Dynamics of strongly driven open systems requires spe- 
cial attention, as the external field modifies not only the 
Hamiltonian of the investigated system, but also its in- 
teraction with the environment. Here we consider a pe- 
riodically driven two-level system that exhibits a series 
of avoided level crossings. This "periodic version" of the 
Landau-Zener-Stiickelberg model can be described by a 
Hamiltonian 

TT r.\ ^ ( acos(flt) 5/2 \ 

Hs{t) = acosma. + ^ 5/2 -aco^m ) ' 

(1) 

where ^l = 2'k/T is the angular frequency of the exter- 
nal field and the second equality holds in the {|+), |— )} 
eigenbasis of ct^. Although similar model Hamiltonians 
can describe a large variety of physical systems, for the 
sake of definiteness we will consider a spin 1/2 system 
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being coupled via Zeeman interaction to an oscillating 
magnetic field. We assume that the spin is embedded in 
a bosonic reservoir (which, in this aspect, can represent 
phonons if the distinguished spin is coupled to a crystal 
lattice) described by the Hamiltonain 



E 



(2) 



where Ofc and a^^ are the annihilation and creation oper- 
ators of the k-ih mode satisfying [afc/,afct] = 5kk' ■ (Note 
that we set h = 1.) The system-environment interaction 
Hamiltonian is written as 



V ^S<. 



' y^gfc(afc^ 

k 



Ofc), 



(3) 



where S can represent any (hermitian) spin operator, and 
the coupling constants gu are assumed to be real. (Note 
that this kind of interaction resembles the coupling of 
a two-level atom to electromagnetic field modes [1^, as 
well as to spin-phonon coupling in solids j29|.) 

Considering only the periodic Hs as the generator of 
the time evolution (no environmental influence) , Floquet 
theory tells us that it is possible to find a time-dependent 
eigenbasis 



(wi(i)|u2(i)) =0, 



\Ur{t + T)) = \Ur{t)), 
{Ur{t)\Ur{t)) = 1. 



(4) 



Unlike the Floquet states \ur{t)), the elements of this 
basis themselves are are not T-periodic functions, which 
is related to a nontrivial phase effect [13]. Let us re- 
call that if is a Floquet quasi-energy and the corre- 
sponding state is |(/)r(t)), then the same holds for tr + n^ 
and \<i>r{t)) X exp(ml7t) for any integer n. However, these 
states are equivalent from the dynamical point of view, 
thus it is sufficient to focus on the two nonequivalent 
quasi-energies the magnitude of which are the closest 
to zero. For the sake of definiteness we will assume 
ei < 62. Note that the Schrodinger equation induced by 
the Hamiltonian ([T]) can be rewritten as an inhomoge- 
neous differential equation of Mathieu-type, numerical 
methods for computing ei and £2 using this fact can 
be found e.g. in Ref. [31|. Having obtained the quasi- 
energies and the corresponding states given by Eq. ([4]), 
the time evolution operator U{t) — \4'r{t)) {4'r{^)\ can 
be constructed. Then, returning to the open system, we 
can transform V into an interaction picture with and 
Hs being the Hamiltonians of the uncoupled total sys- 
tem. In this way, we obtain: 



Si{t) = J2 S{ij)e-"^' + h.c. 



UJ>0 



where 



S{u;)= J2 \urm{u.AOMr'\S\r))r., 



T 



Ur-i {t)\S\Ur(t))dt, 



(5) 

(6) 
(7) 



with n being integer, and the sum runs over indices satis- 
fying e^/ — — f^n = uj [s^ . This means that ~ unless the 
Fourier component S{u!) is zero - three sets of positive 
frequencies appear. 
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nVL. 



(8) 



Assuming that initially the density operator of the com- 
plete system factorizes, p(0) = Ps(0) ® Pr(0), standard 
methods (Born-Markov approximation) lead to the in- 
teraction picture master equation 



^p,(t) = -Tr,y [V{t),[V{t - s), ps{t) ® pr]]ds, 



(9) 



where Tr^ means trace over the reservoir degrees of free- 
dom, and the explicit indication of the interaction pic- 
ture has been omitted. Next we insert Eq. ([5]) with the 
standard periodic (exp(±iWfci)) time dependence of the 
interaction picture creation and annihilation operators of 
the bath into the master equation above. Then the sums 
over environmental modes appearing in Eqs. ^ and ^ 
are transformed into an integral over the mode frequen- 
cies with the mode density ^{bj) being a weight factor. 
Finally we perform rotating wave approximation (RWA) 
in the system-environment interaction (3^ . (Note that 
we did not use RWA for the coupling of the system and 
the external field.) With this approximation the integral 
in Eq. ([9]) can be evaluated, leading to an interaction pic- 
ture Born-Markov master equation for the reduced den- 
sity operator of the two- level system [s^ : 



dps 
dt 



u>0 



X (^S\u)p,S{u) - ]^S{u)S\u:)ps ~ ^PsS{u)S\uj) 



+7'(a;) ( S(u;)psS\u;) - -S\lo)S[lo)ps 



-PsS\lu)S{uj)\ , 



where 



7(w) 



V{uj)g\u:){n{uj)), 



(10) 



(11) 



with {n{Lo)) = Tr^ {a^ {uj)a{ix!)pr) representing the aver- 
age number of excitations in the environmental mode la- 
beled by the frequency uj. Returning to the Schrodinger 
picture. Eq. (|10p turns out to be a Lindblad-type equa- 
tion [34] , but as a consequence of the strong driving field, 
the Lindblad operators have explicit time dependence. 
We note that decoherence effects can also be described 
by the aid of the Feynman- Vernon influence functional 
method [s^ , which has already been applied successfully 
[H, [13] to bistable quantum systems (like the one con- 
sidered here). 

Additionally, as a result of the system-environment 
RWA, the frequencies appearing in Eqs. ^ and (fTU|) 
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are exactly the same. Recalling that if is a Floquet 
quasi-energy, then the same holds for + nfl for any 
integer n, it is possible to establish a connection be- 
tween the Floquet spectrum and the combined energy 
levels of a two-level system and a quantized single mode 
field 1^32]. In view of this, the master equation above 
can be interpreted as the coupling of those transitions of 
the combined spin-field systems to the resonant reservoir 
mode, where the matrix element of the coupling opera- 
tor is nonzero. Technically, it is convenient to collect the 
terms in Eq. (jlOp that contain the same operator part (eg. 
|ui(0)) (1*2(0)1) then to calculate their common coefficient 
by evaluating the sum over frequencies for a sufficiently 
large, but practically finite number of modes (see sub- 
sec. IIII Bp . Note that interaction with the environment 
usually also induces a Lamb-type renormalization of the 
energy spectrum of the system, but in this framework the 
renormalization Hamiltonian (that should appear on the 
rhs. of Eq. (jlOp in a commutator with ps) commutes with 
the system Hamiltonian Hg [32|, thus it has no special 
importance from our point of view. 

Finally, let us note that in the optical case (when the 
offdiagonal elements of Hg oscillate) with system-field 
RWA and exact resonance the Floquet quasi-energies and 
states can be calculated analytically, see Refs. [H, Issj , 
where the quantum jumps leading to the strong-driving 
MoUow spectrum are also introduced. 



III. DIRECTION AND CHARACTERISTIC 
TIME OF DECOHERENCE 

A. Summary of the free time evolution 

Before investigating how the environment induced de- 
coherence modifies the dynamics of the spin system, it is 
certainly worth recalling the main features of its free (uni- 
tary) time evolution. First, let us note that by exchang- 
ing (7x and az in the Hamiltonian ([1]), we obtain a hermi- 
tian operator Hg that is unitarily equivalent to the origi- 
nal system Hamiltonian: Hg = 1/2(0"^ + (Jx)Hg{(Jz -f (Jx)- 
As Hg can describe a two-level atom subjected to elec- 
tromagnetic field (classical Rabi problem without RWA 
[11]), several studies have been devoted to the time evolu- 
tion induced by this Hamiltonian. Floquet analysis of the 
dynamics has already been done in Ref. [38] , while more 
recent results (see [13, [H, HI] and references therein) 
usually focus on the applications in quantum information 
processing. Considering a different field where the model 
above appears, let us note that as an additional effect un- 
derlying the analogy between atoms and quantum dots, 
exciton Rabi oscillations have also been detected [40, 4l[. 
In these strongly confined solid state systems, however, 
memory effects are important, the dynamics can often be 
non-Mar kovian [i^ . 

Depending on the parameters, Hg can induce diverse 
dynamical behaviors, which are qualitatively different. 




FIG. 1: a) The difference of the Floquet quasi-energies je2 — ei | 
for the dimensionless Hamiltoninan given by Eq. (I12|l as a 
function of A and A. b) The LZS parameter Plzs = 1 — 
exp(— 7rA^/2A) as function of A and A. 

In order to reduce the number of the parameters, we 
introduce a dimensionless time variable t = f2t, and ob- 
tain that without environmental effects, the Schrodinger 
equation 

^^1*) = [Acos{T)a, + |*> (12) 

governs the dynamics, with A ~ a/Q, and A = To 
compare with the LZS model, the parameter that charac- 
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FIG. 2: Decoherence-free time evolution of the expecta- 
tion value (iTz) (without decoherence) for different parame- 
ter values. The initial state is the lower eigenstate of a^, 
Ps{0) = \-){-\. 



terizes the dynamics can be Plzs = 1 — exp(— 7rA^/2A), 
which would be the transition probability if the external 
field were linear with a sweep rate equal to the maximal 
one (at the crossing), i.e., A. From a different point of 
view, strong and weak driving also results in qualitatively 
different dynamical behavior, here A/ IS. is the relevant 
parameter: when A is small compared to A, the popula- 
tions in the eigenbasis of Ux exhibits oscillations with the 
(dimensionless) Rabi frequency VIr = A^ + [1 — A)^, 
while a strong driving field A ^ A forces 27r-periodic 
dynamics. Additionally, when A < A, it is clear from 
the optical analogy, that RWA in the system-field in- 
teraction accurately describes the dynamics. Note that 
in our case this approximation (dropping fast oscillating 

counter-rotating terms in Hg, i.e., replacing terms pro- 
portional to COST in the offdiagonals with l/2exp±ir) 
is equivalent to assuming a rotating external field [43| in 
the y — z plane: 



^RWA ^ A/2cos{T)cr, + A/2sm{T)ay. 



(13) 



Fig. [T] shows the difference of the Floquct quasi- 
energies |e2 — ei| for the dimensionless Hamiltonian ap- 
pearing in Eq. (fT^ as a function of A and A. For Hf^^, 
the difference ei — £2 can be calculated analytically and, 
clearly, it is equal to the Rabi frequency flji. As expected, 
RWA is an accurate approximation for ^ ^ A: Around 
driving field amplitudes A = A, the relative difference 
between the quasi-energies with and without RWA is of 
the order of 10%, and it is increasing for larger values of 
A. Note that parameters when e2 — ei = n, n = 0, 1, . . . 
are special in the sense that in these cases the dynam- 
ics has the same periodicity as the driving field. Fig. [1] 
also shows the LZS parameter PlzSi indicating the rela- 



tion between the LZS-type classification of the parameter 
space and the cases of strong and weak driving. Figure 
[5] shows the time evolution of (cr^) for four qualitatively 
different situations, strong and weak driving with both 
A and A being either larger or smaller than the driving 
frequency $1 = 1. The different periodicity and the role 
of small and large LZS parameter can clearly be seen in 
these plots. 



B. Quasistationary solutions with decoherence 

In this section we investigate the steady state solu- 
tions of the master equation pO|) . i.e., the case when 
the time derivatives of the interaction picture matrix el- 
ements are zero. In fact, these solutions do depend on 
time, the matrix elements will not be constants in the 
Schrodinger picture, but this kind of time dependence is 
well known, and, additionally, considering a Hamiltonian 
with explicit time dependence, it is clear that generally 
there are no constant solutions. In fact, steady state so- 
lutions of Eq. (fTO|) in the sense above are periodic (il ] 
due to the time evolution of the Floquet states \ur{t)). 

Now we specify the terms T'luj) and g{Lu) in 7 and 7' 
(Eq. pTjl ) that describe the mode density of the reservoir 
and the strength of the coupling of the spin system to 
the environmental mode with frequency lu. Assuming a 
thermal bath, the average number of excitations {n{Lu)) 
is given by the Bose-Einstein distribution. Recalling the 
case of a two-level atom in thermal electromagnetic field 
[2^ and a molecular nanomagnet in phonon bath [29l |. 
the choice 



7(tj) = 



7 (w) = 



QkTr _ 1 



(14) 
(15) 



is rather general (recall that H = I). Here the oj inde- 
pendent K describes the overall strength of the coupling, 
and Tr is the temperature of the reservoir. Addition- 
ally, it is instructive to rearrange the terms in the master 
equation pU]) by collecting the coefficients of the interac- 
tion picture matrix elements pij = {ui{0)\dps/dt\uj{0)). 
Focusing on the change of the populations, we obtain a 
Pauli-type equation 

Pii = K [P22 (r2i + r'12) - pii (ri2 + Fai)] , 

P22 = "Pii, (16) 

while the offdiagonal elements decay according to 



P12 = P21 



(ri 



ri2 + r2i — 



-rn + rf,2 + r;2 + r^i-2r^)pi2. (i7) 



The coefficients on the rhs. of Eqs. (fTB|) and (|17p are 
determined by the type of the system-environment cou- 
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FIG. 3: The quasistationary expectation value {o-z)'° = 
Tts (p'^CTz) versus the driving field amplitude for A = 1.5 
and different temperatures. 



pling S and parameters A and A in the following way: 

n>0 
n>0 

^ 7(^,7) l((l|^|2))„|', 

Ji>0 

^7(^°)((2|5|2))„((1|5|1))_„, (18) 



T21 



n>0 



and the primed quantities (that vanish at zero tempera- 
ture) can be obtained by the substitution 7^7'. Setting 
the Ihs. of Eqs. pBl[T7|) to zero, the solution is a diagonal 
density matrix (in the |ui(0)), |u2(0)) basis), where the 
ratio of the populations is given by 



pH _ r2i+n, 

OS 

P22 



ri2 + r 



(19) 
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Clearly, p\\+P22 = 1: and the time evolution converges to 
the solution above in the long time limit: pf^ = pii{oo). 
Let us recall that it is possible to interpret + V'21 
(r2i + ^12) the sum of the environment induced 
\u2) (IM2) transition rates over differ- 

ent number of excitations in the driving field. This im- 
plies that for strong driving, the steady state solution 
at zero temperature will not necessarily be the ground 
state \ui) : the coupled system of the spin and the driv- 
ing field can emit excitations into the reservoir by both 
transitions |u2), \u2) provided the net 

energy flow is directed towards the environment. This 
process is possible in the \ui) — > \u2) (from ground to 
excited state) transition as well, since the energy gain of 
the spin system can be compensated by an appropriate 
loss in the energy of the driving field. 



This effect combined with the parameter dependence 
of the quasistationary states can be used for state prepa- 
ration: As an example, Fig. [3] shows the quasistation- 
ary expectation value {uzY^ = Tr^ (p'^'^a^) as a func- 
tion of the driving field amplitude for different temper- 
atures. As we can see, there are several points, where 
the long time limit solution is basically a certain eigen- 
state of CTz at time instants r = 2m:. (For population 
transfer at periodic crossings using a different method see 
Ref. [ii].) Note that the vahdity of RWA in the system- 
environment interaction (this assumption led us to the 
master equation (jlOp . or, in other words, it resulted in 
the fact that there is only a single uj in each term of 
the sums defining the coefficients F in Eq. ((T8|)) requires 
\er - er'\ > 7(u;„) |((i|5'|j))„|^ . Clearly, when ci = €2, 
this condition cannot be met, but the parameter values 
where the quasistationary sates are the eigenstates of CTz 
are far from these degenerate points. 

On the other hand, for weak driving, only the ground 
state \ui) will be populated in the long time limit. In 
this case, apart from resonance A = 1/2, \ui) is basically 
the equal weight antisymmetric superposition of the CTz 
eigenstates, |wi)(0) ~ (|+) — |— ))/\/2, implying that in 
the quasistationary case (cr^) « 0. 

For high temperatures ~ F21, thus, as it is ex- 
pected, in this case the reduced density operator of the 
spin system will always be proportional to unity in the 
long time limit. Let us note, however, that temperatures 
in the mK range can already be termed as high in the 
context above: if H. has the order of magnitude of MHz, 
kTr/fl = 100 is satisfied with « 1 mK. 



C. Dynamics: Decoherence time and pointer states 

The results of the previous subsection regarding the 
quasistationary solutions do not depend qualitatively on 
the type of the system-environment coupling, but for 
dynamical calculations, we have to specify the opera- 
tor S. In the present paper we investigate the cases 

S — CTy, Uz- 

Figure |4] shows examples of dynamics of the expecta- 
tion value ((Tz) for different system-environment coupling 
strengths, temperatures and a coupling operators. As we 
can see, the oscillations seen in the free time evolution 
are damped in this case. The higher the temperature, the 
stronger this damping effect is. As temperature modifies 
the final (quasistationary) state according Eq. the 
amplitude of the long time limit oscillations of {(7^) is 
also temperature dependent. 

In our case Eq. (fT7)) shows that the final reduced den- 
sity operator in the relevant interaction picture will be 
diagonal in the Floquet basis for any type of system- 
environment coupling. In other words, decoherence 
drives the system into an incoherent sum of these states, 
thus they can be considered as pointer states [1^ in our 
system. However, there are important differences from 
the usual picture of decoherence. The first, and prob- 
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FIG. 4; Time evolution of the expectation value (cr^) for 
different parameters and decoherence rates [A = 0.1, A = 
0.5, K = 1.0 for the upper, and A = 10, A = 0.4, k = 10"* 
for the lower graph), system-environment coupling operators 
and temperatures. The initial state is the lower eigenstate of 
Ps(0) = | — )( — |. Note the emergence of solutions with the 
same periodicity as that of the driving field. 



ably less notable point is the fact that in our case it is 
difficult to clearly distinguish between the energy transfer 
between the investigated system and its environment and 
decoherence (which, in this context refers to the bare loss 
of quantum coherence). The time scale of these two, con- 
ceptually different processes is roughly the same (math- 
ematically this is reflected by the fact that the rate of 
change of the diagonal and offdiagonal elements of ps is 
comparable), thus one can not conclude that first fast 
decoherence takes place, which is followed by a slow dis- 
sipative process leading to thermal equilibrium with the 
environment. In some sense it is a size effect: in larger 
systems with more degrees of freedom it is possible to 
make a dynamical distinction between decoherence and 
dissipation based on the time scales (see e.g. Ref. (46|'). 
On the other hand, the nature of the problem implies 
that we have time dependent pointer states: even in the 
long time limit, when ps is diagonal in the Floquet states, 
{az) oscillates as a consequence of the time dependence 
of \ui) and |u2). 

The characteristic time of the decoherence, r^, accord- 
ing to the considerations above, can be defined as the 
time instant when the offdiagonal elements of the reduced 
density operator become smaller then an appropriately 
chosen percentage of their initial magnitude. Recalling 
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FIG. 5: The rate of the decoherence, t^^ as a function of 
the driving field amplitude for different system-environment 
coupling operators and temperatures. 



Eq. pT]) . this definition can be reformulated as 
Td = Re- (Fn + F22 -I- ri2 + T^i - 2T3 



Til + T^22 + ^12 + r21 -^Fg) 



(20) 



Note that the expression above is independent of the ini- 
tial state of the system, it is valid also in the case when 
the offdiagonal elements of the density operator are zero 
at T = 0. Fig. [5] shows tJ^ as a function of the driv- 
ing field amplitude for different system-environment cou- 
plings S. As we can see, the overall tendency is the accel- 
eration of decoherence as the amplitude increases, which 
is related to the width of the distributions ((r'|S'|7-))„, 
i.e., by increasing A, there will be more terms in the co- 
efficients F that are not negligible. The fine structure of 
the curves is determined by the amplitude dependence of 
the Floquet quasi-energies and states, the maxima corre- 
spond to the cases of ei ~ £2. We can also observe that 
different coupling operators induce different decoherence 
rates even for the same value of k. Note that this effect 
can already be seen in Fig. [H where the damping of the 
coherent oscillations were different for S = (Jxi'^yi'^z- 
It is particularly interesting that for the case of A ap- 
proaching zero, Sy and Sz lead to a finite decoherence 
time, while for Sx decoherence effects become negligible 
in this limit. This behavior can easily be explained by 
using the approximate RWA solutions for and |u2) 
that contain basically two Fourier components, and lead 
to vanishing (finite) F coefhcients for Sx {Sy and Sz)- 

In an experimental situation, where {az) can be mea- 
sured with a high enough temporal resolution, and the 
coupling of the system to the environment is unknown, 
the facts discussed above can be used to gain information 
concerning the nature of this interaction by varying the 
amplitude and orientation of the external field. 



7 




F=A cos X 



FIG. 6: Quasistationary magnetization curves for different 
parameter values (A — 1.5, A — 20 for the upper, A = 
0.6, A — 0.1 for the lower curve) and temperatures. Note 
that for small values of A, the Floquet states have only a few 
Fourier components that are not negligible, but this number 
increases for stronger driving. 



IV. HYSTERESIS CURVES: FROM A SERIES 
OF LOOPS TO LADDERS 

It is often instructive to investigate the response of a 
physical system subjected to a periodic field as a func- 
tion of this field itself. Assuming that our two-level sys- 
tem represents a spin 1/2 particle, (cfz) is proportional to 
its magnetization in the z direction, i.e., parallel to the 
dimensionless driving field F — A cost. Based on this 
aspect of the model, and the fact that - as we shall see 
- the functions {(Tz){F) are usually multi- valued, these 
plots can be called hysteresis curves. 

Let us start with the quasistationary solutions of the 
master equation (jlOp . when the density operator is diag- 
onal in the Floquet (pointer) basis, with the populations 
given by Eq. p^ . In this case the periodicity of the Flo- 
quet states imply that the hysteresis curves are closed 
lines. The oscillations of the curves shown in Fig. [H] are 
related to the time evolution of the Floquet states, the 
more Fourier components these pointer states have, the 
more local maxima and minima can be observed in the 
graph {crz){F). 

Since for high temperatures the quasistationary solu- 
tion is proportional to the unit matrix (implying (az) ~ 
0), we expect the functions {az){F) to be squeezed in 
the vertical direction as the temperature increases. This 
effect can clearly be seen in Fig. [HI for high enough tem- 
peratures the plotted curves are close to a horizontal line, 
independently from the parameters. 

Unless the initial density operator of the spin sys- 
tem is equal to the quasistationary solution, the func- 
tion ((Tz) {F{t)) obtained during the whole time evolu- 
tion r = . . . oo is multi-valued, and this hysteresis curve 



FIG. 7: Convergence towards the quasistationary magnetiza- 
tion curves at high temperatures {kTr/il = 600). The grey 
dot and the arrow indicate the starting point (correspond- 
ing to Ps(0) — j— )(— I) and the initial direction. The right 
panel shows the case of weak decoherence, while environment 
induced effects are moderate for the left panel. (Note that 
according to Fig. [H the same value of k means faster deco- 
herence for larger driving field amplitudes.) The insets show 
(az) as a function of t/2tt. 



reflects the convergence of the system to the long time 
limit solution. That is, as Figs. [7] and [H] show, a quasista- 
tionary magnetization curve acts as an attractor, all the 
paths (cTz) {F{t)) converge to this curve independently 
from the initial point (cr^) (-F'(O)) . Visually, the length 
of the path till the quasistationary curve is reached is 
related to the dynamics: if there are a lot of detours be- 
fore getting close to the final curve, decoherence is slow; 
on the other hand, when only a single line is visible to- 
wards the quasistationary curve, coherent oscillations are 
damped strongly. 

Let us note, that in certain experimental situations 
the dimensionless amplitude A can fall orders of magni- 
tude beyond the applicability of the method described 
so far. The main problem is the determination of the 
Floquet eigenstates, which requires the diagonalization 
of a matrix, the dimension of which is proportional to A. 
Therefore we developed an alternative, approximate way 
of calculating the parameters appearing in the dynamical 
equations (fTHl fT7| : The eigenvalue equation that leads to 
the Fourier components 

{±\u,)n= r e^'^'''\±\u,{t))^dt (21) 

J 

can be rewritten as a system of four differential equations 
with n being considered as a continuous variable. These 
equations are coupled via a term which is proportional to 
A/A, which can usually be considered as a small number, 
allowing for the equations to be solved iteratively. The 
zero order approximation (assuming A/ A = 0) is given 
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F=A cos T 



FIG. 8: Convergence towards the quasistationary magnetiza- 
tion curves at low temperatures {kTr/il — 1). The grey dot 
and the arrow indicate the starting point (corresponding to 
Ps(0) = | — )( — I) and the initial direction. The insets show 
{(Tz) as a function of T/2n. 

in terms of the Airy functions 

(-K)„ = iV_Ai(^(^|y\n-A)^ , 

(+|U2)„ = {-\ui)n, {-\u2)n = -(+|mi)„, (22) 

where N± are constants. The next iteration provides so- 
lutions where the oscillatory part of the Airy functions 
(for negative arguments) are also strongly damped for 
|n| > A. Using these analytical approximations, we found 
that in the high temperature limit, the coefficients F are 
proportional to A^ (indications to this kind of behavior 
can already be seen in Fig. [5]) . Note that these analyt- 
ical results were verified by numerical calculations with 
the largest amplitudes allowed by our computational re- 
sources. In fact, not only the scaling with were seen, 
but a good agreement concerning the prefactor has also 
been found. Similarly, the approximate eigenstates (|22)) 
in those domains where the Airy functions do not oscil- 
late too fast were also found to be in a reasonable agree- 
ment with the numerically exact results. 



Application to molecular nanomagnets in slowly 
oscillating, large amplitude external magnetic fields 

At high temperatures, when (cr^) = is the final solu- 
tion, the hysteresis curves can reach the final horizontal 



line via a series of steps (see the left panel of Fig. [7]). 
According to a recent experimental result [43] , this kind 
of behavior can be relevant in physical systems where 
crystals consisting of high-spin molecules such as Mni2- 
Ac (or simply Mni2) and Fcg (also known as molecular 
nanomagnets are being driven by periodic external 
magnetic fields. These special molecules contain transi- 
tion metal atoms with strongly exchange-coupled spins, 
which causes them to behave as a single, large spin. Ex- 
periments on the magnetization dynamics of these molec- 
ular crystals have shown the presence of a series of steps 
in the magnetization curve at sufficiently low tempera- 
tures [H, l4^, [sO, HH . This behavior is a consequence 
of quantum mechanical tunneling of spin states through 
the anisotropy energy barrier and occurs when the exter- 
nal field brings two levels at different sides of the barrier 
into resonance via Zeeman interaction. When the exter- 
nal magnetic field is swept linearly, an appropriate LZS 
model around a certain resonance provides a very useful 
approximate description of the dynamics (physical conse- 
quences of the difference between a model involving not 
only two energy levels and the LZS treatment can be 
found in Ref. [ISj). 

In the following we focus on the molecule Mni2, which 
can be considered as a representative example of molec- 
ular nanomagnets, and, consequently, it has been inves- 
tigated in several important experimental works, includ- 
ing the one reported in Ref. [47|. In this experiment the 
sweep rate was 5.83 mT/s, with amplitudes around 0.25 
T, while the temperature was = 0.25 K, corresponding 
to the high temperature limit kT^/Ml 3> 1. As an exam- 
ple, let us concentrate on the seventh resonance around 
3.67 T, where the magnetic levels labeled by m = — 10 
and to' = 3 correspond to the cr^ eigenstates |— ) and 
\+) in our model, respectively. The level splitting at this 
anticrossing can be calculated using the appropriate spin 
Hamiltonian [l^H^, leading to 7 x 10~^ K in temper- 
ature units. Combining these values, the dimensionless 
parameters are A « 10^^, A « 10^, thus the approxi- 
mate method described previously in this section has to 
be applied. 

Physically, the overall decoherence rate (resulting from 
various effects [l^) in Mni2 is around 10^-10^ 1/s, which 
means an extremely fast process compared to the oscilla- 
tion of the external field, where S7 « 0.1 1/s. If we assume 
that the phonon bath is the only source of decoherence 
(which is not the case, see below), it is impossible to 
obtain ladder-like hysteresis curves similar to the experi- 
mental results. For very strong phonon induced decoher- 
ence, (ctz) as a function of the external field would show a 
sudden convergence towards its stationary value of zero, 
in other words, the hysteresis curve would consist of two, 
almost perpendicular straight lines, starting with a verti- 
cal one that connects the initial point and the stationary 
horizontal line. (Similarly to the curve corresponding to 
the strongest decoherence in Fig. [51) However, both ex- 
perimental and theoretical results show (see eg. Refs. (52j 
and [1^) that the main source of decoherence is related 
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FIG. 9: Hysteresis curves for slowly oscillating external fields 
calculated with assuming incoherent transitions as a conse- 
quence of dephasing. The figure corresponds to the seventh 
resonance (around 3.67 T) of the molecular nanomagnet Mni2 
driven by an external field with amplitude of 0.25 T and 
Q = O.f 1/s. The ladder-like curve for the weakest phonon 
induced effect (dotted line) reflects qualitatively the results of 
Ref. El. 



to dipolar and hyperfine interactions which modify the 
local environment of the spins and eventually result in a 
distribution of the level splittings. Investigating the dy- 
namical properties of our system, it can be seen (both an- 
alytically and numerically) that this effect decreases the 
ofFdiagonal elements of the density matrix in the Floquet 
basis. The characteristic time of the process depends on 
the width of the distribution of the level splittings, but 
if we take into account that a single period of the ex- 
ternal field takes roughly 10 s, dephasing due to dipolar 
and hyperfine interactions is practically instantaneous. 
We note, however, that this kind of dephasing usually 
does not lead to a horizontal stationary magnetization 
curve, as even the incoherent sum of the Floquet states 
has nontrivial time dependence. 

This implies that it is worth assuming that the re- 
sults reported in Ref. [i^l reflect the interplay of fast 
dephasing due to the distribution of the level splittings 
and a much slower phonon induced process. Therefore 
we calculated the dynamics with constantly zero offdiag- 
onal density matrix elements in the Floquet basis (corre- 
sponding to instantaneous dephasing) by the aid of the 
approximations that led us to Eqs. ip^ . As it is shown 
by Fig. [21 the value of = 0.005 1/s leads to qualitative 
agreement with the experimental results. However, it is 
clear that decoherence rates F are not variable parame- 
ters, their values are determined by the physical system 



we are considering. Concretely, the value of k in Eq. [15] 
(that can be calculated according to Ref. [2§]) and the 
scaling of T{A) as together determine the decoher- 
ence rates. In this procedure there are basically no free 
parameters (apart from some uncertainty of the sound 
velocity in the molecular crystal) and we obtained that 
is around 0.025 1/s, which, according to Fig. [51 is 
slightly too large, it leads to hysteresis curves with a few 
steps only. However, keeping in mind the approximate 
nature of our treatment (considering not only the calcula- 
tion of the Floquet states , but also the spin-phonon 
couphng operator, see eg. Ref. [2§]), we think that the 
methods developed in this paper provide a solid starting 
point for a theory that aims to describe nanomagnets 
driven by periodic magnetic fields in a quantitative way. 
Concerning our current results, we can conclude that the 
thermal phonon bath alone can not be responsible for 
the experimentally observed ladder-like hysteresis curves, 
but taking the strong dephasing also into account, this 
kind of behavior can be explained. 



V. CONCLUSIONS 

We investigated a two-level system which is driven by 
periodic external field and which is also in interaction 
with a thermal bath. It has been found that - indepen- 
dently from the decoherence rate, from the initial state 
and even from the type of the system-environment cou- 
pling - the time evolution is directed towards an appro- 
priate incoherent sum of periodic Floquet states. The 
final ratio of the populations related to these time depen- 
dent pointer states is determined by the parameters of the 
system Hamiltonian, the type of the system-environment 
coupling and the temperature. Our results show that 
the form of the quasistationary hysteresis curves is com- 
pletely determined by the time evolution of the Floquet 
states, while the rate of convergence towards these curves 
is related to the characteristic time of the decoherence. 
As an important example of the possible applications, 
we have shown that our model can be used to describe 
molecular nanomagnets driven by periodic external mag- 
netic fields. 
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